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We explain excess events near the endpoints of the double beta decay {(3(3) spectra 
of several elements, using the neutrinoless emission of massless Goldstone bosons. Models 
with scalars carrying lepton number —2 are proposed for this purpose so that ordinary 
neutrinoless is forbidden, and we can raise the scale of global symmetry breaking above 
the 10 keV scale needed for observable emission of conventional Majorons in (3p. The 
electron spectrum has a different shape, and the rate depends on different nuclear matrix 
elements, than for the emission of ordinary Majorons. 
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One of the deep questions of particle physics is whether there exist any fundamental 
particles of spin zero. A candidate is the Majoron, the massless Goldstone boson that 
would exist if lepton number (L) were spontaneously broken . Majorons would be hard 
to detect since they couple directly only to neutrinos, with a strength proportional to the 
neutrino masses divided by the L-breaking scale v: 



If Qi, were large enough, Majorons could be seen in double beta decay a rare 

process now observed in seven elements. Besides the usual neutrino-emitting mode (3^2^, 
a Majoron could be emitted through annihilation of the virtual neutrinos, (3 (3m- The (3 (3m 
and (3(32u signals are distinguishable since for (3(3m the decay energy is shared among fewer 
particles, skewing its electron spectrum toward higher energies. 

In fact, a mysterious excess of high-energy electrons is seen in the (3(3 spectra for 
several elements. Such an observation was first made in 1987 for ^^Ge ^^^Se +2e- by 
Avignone et al. [|| , although they, and other groups, subsequently excluded a signal at the 
original level |]5| . Now the UC Irvine group also finds excess numbers of electrons near but 
below the endpoints for ^'^'^Mo, ^^Se and ^^''Nd, with a statistical significance of 5(t 0. 
Such events also persist in ^^Ge [0, at approximately a tenth of the original rate. 

Formerly all these (3(3 measurements were compared with the Gelmini-Roncadelli (GR) 
model H] 0], which however has been ruled out by LEP's bounds on the invisible Z width. 
In this model lepton number was spontaneously broken by an electroweak-triplet Higgs 
field, resulting in both a Majorana mass and a direct coupling to the Majoron for u^,, 
as in eq. (j^). The presently observed excess events could then be explained by (3 (3m if 
g,^^ ~ 1 X 10~*, and the absence of (3(3qi, decay — which would appear as a line at the 
endpoint in the sum-energy spectrum — requires ttIj,^ < 1 eV. Eq. (|lD then implies that 
the triplet VEV must be unnaturally small. 



(A similar bound follows from astrophysical considerations.) 

With the demise of the GR model, it appeared that there existed no models capable 
of predicting (3(3m at an observable rate. In this Letter we propose an alternative broad 
class of Majoron models which might explain the excess events, while still preserving the 



9u = m^/v. 



(1) 



^; < 10 keV. 



(2) 
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agreement between theory and geochemical experiments for the decay rate of ^^^Te, 
i30Te and ^ssu §. 

We start by describing some generic features of PPm decays, and then turn to model- 
specific issues. First, agreement with the Z-width measurement rules out a direct Majoron 
coupling to the Z boson. Lepton number breaking, if it occurs, must come from an 
electroweak-singlet field. Majorons therefore cannot have renormalizable couplings to Uf,; 
instead, (3(3m proceeds through mixing of Uf, with sterile neutrinos that couple to Majorons. 

Next one must explain how (3(3m but not /3/3oiy events are seen, since lepton-breaking 
generically gives both effects. The GR model did so by requiring a small lepton-number 
breaking scale as in eq. (^, and the same expedient is used by alternate models where 
electroweak-singlets rather than triplets do the lepton-breaking. But once such a small scale 
is introduced by hand it is unnecessary for the scalar emitted in /3/3m to be a Goldstone 
boson; it could just as well be given a mass of 10 keV. We refer to these as 'ordinary 
Majorons' (and denote their emission by PPom) since they share many features (electron 
spectra, etc.) of the GR model. These models all suffer from a hierarchy problem since 
the lepton-breaking VEV is not stable under renormalization and so must be artificially 
fine-tuned [1T0 |. 



Here we introduce a class of models that do not require such a small scale, inspired 
by realizing that the smallness of v in eq. (^ comes from the bound on P/Sqi, decay. This 
restriction disappears if a lepton number, L, carried by electrons, remains unbroken, so 
that /3/3oi/ is completely forbidden. Then (3(3m occurs only if the scalar particle itself is 
'charged,' i.e., it carries L{(f) = —2, in units where L{e~) = +1. Fine-tuning of the scalar 
mass is avoided if the scalar is a Goldstone boson, a 'charged Majoron.' All of this can 
be realized in a simple way, using for example a global SU{2)s x U{1)l' symmetry broken 
down to U{1)l- (The U{1)l' factor exists so that Majoron couplings to ordinary charged 
leptons can be avoided.) Two of the resulting Goldstone bosons then carry the unbroken 
charge, just as the would-be Goldstone bosons that make up the longitudinal VF^'s carry 
electric charge in the Standard Model. We call all such models 'charged' Majoron models 
and denote the resulting process (3(3cm- 

In what follows we compare charged and ordinary Majoron models with each another 
and the data. We will: (z) show how the models can be experimentally distinguished by 
the predicted shape of the sum electron spectrum, and (ii) find the necessary conditions 
for getting a large enough rate of /3/5m • 



i. The electron spectrum. The decay rates for /3/92i/, PPom and /3/9cm all have the form 

drm = ^ \A{pp)fdnr., (3) 

where Gp is the Fermi constant, ^ is a matrix element (more about which later), and dfln 
is the phase space, a function of the momenta and energies of the two outgoing electrons, 

1 ^ 
dn^ = ^^{Q - - 62)" Yl Pk€kF{ek)dek. (4) 

k=i 

-F(e) is the Fermi function, equal to 1 when the nuclear charge vanishes, and Q ~ 1 — 
2 MeV is the endpoint energy for the final-state electrons. 

The integer n in eq. @j depends on the decay channel: n = 5 for (3(32^, 3 for (3 (3cm, and 
1 for /3p oM • ^ controls the shape of the electron energy spectrum, because for these decays 
all lepton energies are <^ pp ^ 100 MeV, the nuclear scale that determines the matrix 
element A. Thus it is a good approximation — relative error of order Q/pp ~ 1% — to 
neglect ^'s dependence on the outgoing lepton energies. In particular, near the endpoint 
the electron spectrum in (3(3cm goes like (Q — E)^ as E, the sum of electron energies, 
approaches Q. This is intermediate between {Q — E) for (3f3oMi and {Q — E)^ for (3(32u- 
The spectra are plotted in Fig. 1. The distinction between (3 (3cm and (3(3om is one of our 
main results. 



Figure 1. The electron sum-energy spectrum for (3(32^ (dashed), (3 (3cm (solid) and 
(3(3oM (dot-dashed) decays. 
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We emphasize that the spectral shape difference characterizes these two classes of 
models, independently of any model-specific details. This traces to the vanishing of (3(3cm 
amplitudes at zero Majoron momentum, q = 0, due to the fact that Goldstone bosons are 
always derivatively coupled. The puzzle is why (3(3om amplitudes do not also vanish at 
g = 0. There one finds that graphs with Majoron Brehmsstrahlung off the electron lines 
survive as g — because the internal electron propagator diverges in this limit, leaving a 
finite result. The same cannot happen for (3 (3cm because the electron-Majoron coupling is 
forbidden by the unbroken lepton symmetry in these models. 

a. The decay rate. Computing the absolute rate requires knowing the matrix element, 
A, of eq. @. Consider arbitrary Yukawa couplings between a set of neutrinos and a 
complex scalar field, (p, 

Cm = Vi{aijPL + hijPn) Vj ip* + c.c. (5) 

Goldstone boson couplings are included here as a special case: if / is the decay constant, 
Xij the conserved charge and rriij the left-handed neutrino mass matrix, then a^j takes 
the form a = —i[X'^m -f mX)/ f, and similarly for bij. Real scalar fields are also easily 
incorporated by taking 6 = a* in all results. For example the GR Majoron coupling arising 
from a triplet interaction C = LeP^TaLe A" + c.c. would give 6*^ = a^e = ig- 
To lowest (zeroeth) order in the scalar energy we get eq. with n = 1 and 



(6) 



Vei is the mixing matrix element for finding the mass eigenstate Ui in an electron charged- 
current weak interaction. The form factors Wa parametrize the matrix element of the 
hadronic weak currents between initial and final nuclei and A^', 

Wc^f, = j d^x {N'\T* [Mx)JpiO)] \N) e'P- 

+ tapapU'^p'' + PaPf3, 

and are functions of the invariants and u-p. is the four- velocity of the initial and final 
nuclei, which are equal since we neglect all dependence in WafS on the final-state momenta. 



Although eq. suffices for (3(3om it happens that for (3(3cm this lowest-order result 
vanishes, making it necessary to work to next-higher order in the outgoing Majoron mo- 
mentum, and leading to the difference in electron spectra (above). The result has the form 
of eq. (I), but with n = 3 and = |W(/3/3cm)|' + |V(/?/3cm)|', where 



^ J [271 p [p^ + mf,. - ie){p^ + mf,. - IS) 



d^p p^{w2 — p^w^] 



(8) 



(27r)^ (p2 -f- ml. — ie){p'^ + ml. — ie) 



Notice that whereas A{I3[3om) has the same combination of form factors as /3/3o!/5 
(3(3cM depends on different nuclear matrix elements. (For I3(3q^, A oc Yli^ei'^fi 
/[(i^p/(27r)^] [(4wi — W2 + P^w^) / {p^ + m^. — is)]] these are related to the Fermi- and 
Gamow- Teller-type form factors by wi = ^Wgt and W2 = Wp + \Wgt-) 

A striking feature of both eqs. (^) and (§) is that they vanish if all neutrinos share a 
common mass, for in this limit the sum over states gives 'Yliij VeiVejttij = aee- The latter is 
zero since the width constraint implies that the Majoron cannot directly couple to Ue. 
As an important corollary, the amplitudes must also vanish if all of the Ui are negligibly 
light compared to the Fermi momentum pp, above which the integrals in (^ and (P) are 
effectively cut off. Thus in any model explaining the excess electron events by Majoron 
couplings to a virtual neutrino, at least one of the i^i must have a mass m^^ > 100 MeV. 

To make contact with the literature, we use the predictions of the GR model as a 
benchmark. In this case 

2 /2 

A{/3(3m)gh = ^glge^M, (9) 

TT 

Qes is the GR coupling of (/? to z^e, and Qa are the axial and vector couplings of the 
nucleon weak currents, and M. = ([a^ ■ am — {9v/gA)'^]h{r)'^ (see eq. (18)). We integrate 
the phase space dfli above a threshold i?th where the anomalous events begin and the 
contribution from ordinary decay should be small. Table 1 compares this with the 
excess events seen by the Irvine group for the elements ^^Se, ^'^'^Mo and ^^*^Nd and in 
the published spectrum of ^^Ge W^■ In all of these cases the excess events comprise R = 2 



to 3 % of the total number observed. We also analyze the geochemically observed decays. 



Table 1: parameters for emission of GR Majorons in double beta decay. T^^J^ inverse 
half-life of the anomalous events; note that for the last three elements this is assumed to be the 
entire rate; |A^p is the Ou matrix element extracted from Staudt et al.\ and R is the ratio of 
anomalous to all events — hypothesized to be 1 for the geochemically observed decays. Oim^ is 
the total phase space for Majoron events (see eq. @), while -Eth (MeV) denotes our choice for 
the threshold value of the sum of the electron energies, above which essentially only excess events 
appear. AOim^ is the part of phase space occurring above -Eth ^-nd fi'eff = 9ee is the coupling 
needed to explain the rate using Majoron emission in the GR model. 



Element 




|A^|2m-2 


R 


Oi 


-E'th 




fi'eff 


^^Ge 


2 X 10-23 


1 X 10^ 


0.02 


2.0 


1.5 


0.9 


1 X 10-* 


82Se 


2 X 10-22 


9 X 10* 


0.03 


17 


2.2 


7.6 


8 X 10-5 




3 X 10-21 


2 X 10* 


0.03 


34 


1.9 


22 


4 X 10-* 


150Nd 


3 X 10-20 


1 X 10^ 


0.02 


260 


2.2 


155 


2 X 10-* 


i28Te 


5 X 10-25 


7 X 10* 


1.00 


0.23 


0.0 


0.23 


4 X 10-5 


i30Te 


1 X 10-21 


5 X 10* 


1.00 


30 


0.0 


30 


1 X 10-* 


238-Q 


5 X 10-22 


3 X 10* 


1.00 


33 


0.0 


33 


2 X 10-* 



assuming all the events are due to (313m- Remarkably, g^Q lies in the range 4 x 10-^ to 
4 X 10-* for all seven elements. 

The predictions for Te are of particular interest because of a recent measurement of the 
ratio of decay rates C = r(i30Te) /V{^'^^Te)= (2.41 ± 0.06) x 10^ [[l|]. If the same couphng 
needed to account for the endpoint anomalies ((/eff ~ 1 x 10-*) is used for these decays 
then the GR and /3/3om models predict too small a ratio: C{(3(3om) = (30.4/0.23) (5/7) = 93. 
The same is not true for (3(3cm models because the rates scale with an extra factor of 
relative to those for (3(3om, and ^'^^Te has a small endpoint energy, Q ^ 0.9 MeV. Thus we 
find CiPPcA-i) = 770, in much better agreement with the experimental value. 

To see whether the ordinary or charged Majoron models can predict a large enough 
rate, we display the combinations of couplings, defined in (H), that are equivalent to the 
GR model with coupling g^s, so far as the total rate is concerned. For ordinary Majorons 

^eff = 2^ {^-^-j ( (p2 + ^2J(^2 + ^2 ) ) • (10) 

Here p is the four-momentum of the virtual neutrino, which as indicated by ( ), is integrated 
weighted by the nuclear form factors. An estimate for the size of g^s is therefore obtained 
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by replacing by the characteristic momentum scale, p| ~ (100 MeV)^. For charged 
Majorons, in order of magnitude 



The two qualitatively new features here are: (1) the additional suppression by the endpoint 
energy, Q/p, due to the softer electron spectrum of these models, and (2) a ratio, X, of 
the l3[3cM to /9/9oM nuclear matrix elements. Later we give a more quantitative comparison 
of these matrix elements, in a specific model. 

Consider now a representative model for each class. For (3(3oAn augment the Standard 
Model with two Majorana-Dirac sterile neutrinos, and s_, and a complex singlet scalar, 
(p. The model consists of interactions that preserve electron-type lepton number, assuming 
Le{PLS±) = ±1, and Le{f) = —2. To satisfy the /3/3oi/ constraint we take {(p) = and 
fine-tune the (p mass to be < 1 MeV as required by the /3/3m kinematics. The spectrum 
contains one Dirac neutrino (\1/h, ^^l) with mass M, "i! ^ = si and = sin 6 + s.^. cos 9. 
The massless orthogonal combination u'^ is predominantly z^e- The remaining two neutrinos 
and Uj- acquire Majorana masses at dimension five. 

In terms of the relevant Yukawa couplings to p, 

= -^9+ {s+PlS+) p - ^g- (s-Pr^S-) p* + c.c. (12) 

the effective /3/3m coupling in this model becomes 

_ / 1 \ / (gV-^+M^cos^g)M^sin^g \ 

Since the light scalars couple to the massless neutrino u'^ through mixing, I3(3om may 
proceed via the exchange of z^g, with the result is nonvanishing as M — * oo with 6 
held fixed. (Decoupling is not violated since 6 ~ 1/M if dimensionless couplings in the 
lagrangian are held fixed for M large.) Using (^p^^ ~ (100 MeV)^ and 6 ~ 0.1 we find 
(7efi ~ lO"'*, as suggested by the (3(3 anomalies. We choose M above the meson mass; 
otherwise a heavy neutrino with such large mixing would have been seen as a peak in the 
K ^ ef fl^ or TT — > ez^ [|4| decay spectra. 



For charged Majorons we suppose that the lepton symmetry group surviving at elec- 
troweak energies \s Gl = SU{2)s x U{1)l', and introduce the sterile Majorana neutrinos 
= (A^+, A^-), s+, and s_, whose left-handed components transform respectively as (2, 0), 
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(1, 1), and (1, —1) under Gl. A sterile scalar doublet 4> transforming as (2, 1) is also added, 
whose VEV breaks Gl —>■ U{l)e where Lg = 2Ts + L' . The renormalizable interactions pre- 
dict three massless neutrinos z/^, i/^ and 1/^, and two heavy Dirac neutrinos of mass M and 
A, say. Uf, is mostly u'^, but mixes with A^_|_ and similarly to the previous model. There 
are three massless Goldstone boson states: a complex scalar ip having charge Lf,{f) = —2 
and a real scalar with = 0. 

For P(3cM the unbroken lepton number not only forbids contributions from graphs 
involving purely massless neutrinos, but it also causes the contribution from each heavy- 
neutrino line to go like rather than 1/M when M > 100 MeV. The equivalent GR 
coupling of the Majoron in /3(3cm is roughly 

if M ~ A is the heavy neutrino mass and v is the G^-breaking scale, (cf. eq. (|Tl]), which due 
to the relation between bij and the mass matrix does reproduce (|T^), despite appearances.) 
For M ~ ~ pp, 9 ~ 0.1 and X ~ 1, ges is of order 9'^Q/pp ~ 10~^ as desired. 

To be more quantitative, we must evaluate the new nuclear matrix elements for PPcm, 
which differ from the usual ones because the leptonic part of the amplitude is proportional 
to the Majoron and virtual neutrino momenta through the factor The po piece (V 

of eq. (H)) vanishes in the integral / dpo, up to small 0{Q/p) corrections. To have a 
0"*" 0"*" nuclear transition, the pi piece {U in eq. (^)) must combine with odd-parity 
nuclear operators, which come from the recoil corrections to the nucleon currents and p- 
wave Coulomb corrections to the electron wave functions. Although they are formally of 
higher order in a or f/c of the nucleon than the usual Gamow- Teller and Fermi matrix 
elements, explicit calculations by nuclear theorists show that they need not all be small, 
as we will show [ p!5| . 

For small mixing angles the matrix element becomes 

in which the nuclear matrix elements are 
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\Mi\ 



9 



9 / _ \ 2 
|2 A ' 



{A, + iA2)h) ; 1-^21= - {Ash) , (16) 



and Ai are operators 

Ai=r- [{CmO-n - CnO-m) + {gv/gAf{^n " D, 



aZr r, , ,o ^ -.1 

(17) 

A2 = -Z-^ \ [9v/9Aj + CFm- CFn- ^CFm " T (Jn " r I ; 
A3 = f ■ [{CmCin + Cr^Om) + (fir^ /S'a ) ^ (Dn + Dm)] ■ 

We use the notation (OnOm) = {0~f\ J2n m '^nOmT~^T^\^t) ^ where r+ is the isospin raising 
operator for the nth nucleon. Definitions of C and D can be found in the reviews by Doi 



et al. [0 or Tomoda [jT^. The neutrino potential h depend on the internucleon separation 
vector r = r„ — and the heavy neutrino mass M: 

where = {p^ + M^)^/^ and /U ~ 10 MeV is the average nuclear excitation energy. 

In Table 2 we show the value of the nuclear matrix elements, (^^ p)^/'^, needed for 
charged Majoron emission to account for the anomalous events, assuming that 9'^{M/v) = 
3 X 10~^ (as would be the case for a neutrino with 17% mixing, the limit coming from 
fiavor universality fl^). For comparison, we also show a similar matrix element that has 
been calculated by Muto et al. il^, namely 



Mn = {9v/9a) (pI) ^ ((D„ xam + cTnX D^) ■ f dh/dr) , (19) 

Except for a minor difference in the neutrino potential, the same matrix element is con- 
tained in the Ai term of eq. (^B|). This is only meant to be indicative because there are 
further contributions to our which do not appear in the literature, and which 

are conceivably bigger than Ain.. 

From table 2 we see that charged Majoron emission explains the data remarkably 
well, if the as-yet uncomputed matrix element for ^ss^j jg gi^ilar to those of the other 
nuclei. In particular the f3(3cM rates for the geochemically observed decays of ^^^Te and 
^■^'^Te dominate over the contribution of (3(32^- As was described earlier, this is in better 
agreement with the data than is the result for (3(3om models. 

Because of the extra suppression by {Q/PfY f^f^cA-n we needed a heavy neutrino 
\1/ mixing with Ue at the 9 = 10% level, with a mass near 100 MeV. But there are strong 
limits on 9 for such a heavy neutrino. Of these, beam-dump experiments do not apply 
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Table 2: phase space and matrix elements for charged Majoron emission in double beta decay. 
The phase space factors are defined analogously to those in table 1. We give the values of the 
matrix elements which are needed to explain the data (assuming all events were PPcm for ^^^U 
and Te) for O^M/v = 3 x 10"^, and values of some representative matrix elements that have 



been calculated; see eqs. (16) and (1£), respectively. 



Element 






(EJ-^^I')'^' needed 


\Mh\ (Muto et al.) 


^^Ge 


4.7 


0.43 


1.1 


1.1 


82Se 


78 


7.9 


0.99 


0.95 




160 


42 


1.4 


1.1 




1570 


320 


1.6 


1.3 


i28Te 


0.13 


0.13 


0.32 


0.92 


isorpe 


110 


110 


0.48 


0.78 


238-Q 


35 


35 


0.61 


? 



to our model because \1/ decays invisibly into a light neutrino plus a Majoron. However, 
searches for peaks in the spectra oin, K ^ ue [|2^ rule out a long-lived heavy neutrino with 
such large mixing, unless M > Mk ~ 500 MeV. For such large M the 1313cm rate becomes 
suppressed, unless the unknown nuclear matrix elements turn out to be considerably larger 
than those of ref. [jl^ . 

But there is a loophole in the peak search bounds: they do not apply if the 100 MeV 
neutrino \E' is a broad resonance. A width of 10 MeV is sufficient to hide the tt 



peak beneath backgounds [^, and the K decay searches do not extend below 140 MeV. 
Such large widths occur in our model if the neutrino-Majoron coupling is strong, g ~ 47r, 
in which case we must repeat the (313m analysis with new form factors parameterizing the 
strongly-coupled neutrino sector. Then dimensional analysis again implies that the rate is 
maximal if the scale for new neutrino physics is of order 100 MeV. It would be fascinating 
if double beta decay experiments gave the first hint of a strongly-coupled neutrino sector! 

We note that f ~ M ~ (few hundred MeV) is also the scale at which the global 
symmetry Gl must break. While this is still annoyingly small compared to the weak 
scale, it is a significant improvement on the much smaller 10 keV scales needed in ordinary 
Majoron models. 

A final challenge confronting any model is the strong limit, Ni, < 3.3, p2[| on the num- 
ber of neutrino species at nucleosynthesis. The large coupling implied by the PP anomalies 
means that Majorons had an undiluted thermal energy density at MeV temperatures, a 
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potential disaster, since each Majoron counts as 4/7 of a neutrino. The three Majoron 
states of the /3/3cm models therefore preclude more than two light neutrinos at nucleosyn- 
thesis, which could happen if (or f^) decays before 0.8 MeV. One possibility is that 
rriv^ ~ 25 MeV (large enough to allow it to decouple before nucleosynthesis but still below 
the laboratory mass bound), with a lifetime r < 10'^ sec |j2^. But even if is light, it can 
effectively contribute as negative number of neutrino species if it decays into with life- 
time 6 X 10~^ sec < r < 2 X 10~^ sec [^, an even better prospect for our models. The first 
scenario, for example, occurs in our ordinary Majoron model via '^/i + mediated 

by a dimension-seven operator {L^H)P^[LT-H)ip*(p. Similarly, we can use dimension-five 
and dimension-seven interactions to give majorana masses and lifetimes to and Vr in 
the charged Majoron model. 

In summary, we have presented two classes of models for Majoron emission in double 
beta decay which may be able to explain excess events near the endpoints of several 
elements without seriously spoiling the agreement of geochemical observations with 
theory. We propose a new class of models involving a 'charged' Majoron that carries 
an unbroken lepton number, which predict a softer sum energy spectrum than that of 
ordinary Majorons, and thus make them experimentally distinguishable from the latter. 
Charged-Majoron models are highly constrained, and require a larger, 100 MeV scale of 
symmetry breaking than the unnaturally small 10 keV scale of ordinary Majorons — a 10^- 
fold improvement. The detailed predictions for these models depend in part on the size of 
certain nuclear matrix elements (eq. (|16|)) which differ from those appearing in the usual 
amplitudes for (3132^, f^f^oi,, or (3(3om- All our models point toward a sterile neutrino with 
a mass of a few hundred MeV, a large (~ 0.1) angle for mixing with Ue, and possibly a 
large (~ 10 MeV) width. Nucleosynthesis requirements suggest that the neutrino which is 
predominantly Uj- should either have a mass that is not much below the laboratory lower 
limit, or else a lifetime < 10"'^ sec. 

We warmly thank F. Avignone, M. Moe and A. Turkevich for information about their 
experiments, T. Kotani, R. Shrock and E. Takasugi for helpful correspondence, W. Haxton, 
M. Luty, N. de Takacsy and P. Vogel for helpful discussions and R. Fernholz for his con- 
tributions, including the figure. This work was supported in part by the Natural Sciences 
and Engineering Research Council of Canada and les Ponds P.C.A.R. du Quebec. 
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